Interaction model for magnetic holes in a ferrofluid layer 
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Nonmagnetic spheres confined in a ferrofiuid layer (magnetic holes) present dipolar interactions 
when an external magnetic field is exerted. The interaction potential of a microsphere pair is 
derived analytically, with a precise care for the boundary conditions along the glass plates confining 
the system. Considering external fields consisting of a constant normal component and a high 
frequency rotating in-plane component, this interaction potential is averaged over time to exhibit the 
average interparticular forces acting when the imposed frequency exceeds the inverse of the viscous 
relaxation time of the system. The existence of an equilibrium configuration without contact between 
the particles is demonstrated for a whole range of exciting fields, and the equilibrium separation 
distance depending on the structure of the external field is established. The stability of the system 
under out-of-plane buckling is also studied. The dynamics of such a particle pair is simulated and 
validated by experiments. 

PACS numbers: 75.50.Mm, 82.70.Dd, 75.10.-b, 83.10.Pp 




I. INTRODUCTION 

The dynamic properties of so-called magnetic holes in 
ferrofluid layers has been the obiect of increasinginterest 

aliiiiii 
12^ . These 

systems consist of spherical non-magnetic particles in a 
carrier ferrofluid, whose size is order of magnitudes (1-100 
jim) above the one of the magnetic particles (O.Ol^m) in 
suspension. The ferrofluid appears then as homogeneous 
at the scale of the large particles - holes -, and their 
effect on the magnetic fleld can be modelled as a dipo- 
lar perturbation, where the magnetic moment is oppo- 
site to the one of the displaced ferrofluid Q . The system 
is generally conflned between glass plates in quasi-two- 
dimensional layers, whose thickness slightly exceeds the 
diameter of the holes. The induced dipolar interactions 
give rise to a rich zoology of physical phenomena, such as 
crystallization of magnetic holes in constant or oscillating 
magnetic fields fEQJ10,0|, order-disorder transitions in 
those crystals or non-linear phenomena in the 

dynamics of those systems in low frequency oscillating 
fields finL Inl Il2. uM ■ commonly described using braid 
theory ^illlulEi El ■ 

The understanding of these systems is important in re- 
lation to industrial ferrofiuid applications [24, 25J, or for 
their potential use in biomedicine [211^12^. The dy- 
namics of these phenomena can also be used indir ectl-y 
to characterize the ferrofiuid's transport properties l2£ 
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as its viscosity. Eventually, the ability to shape the ef- 
fective pair interaction potentials through the imposed 
external magnetic field, makes these system good candi- 
dates as large analog models to study phase transitions 
aggregation phenomena [19] or fracture phenomena 
in coupled granular/fiuid systems. 

Nonetheless, despite the theoretical studies on similar 
dual systems as ferromagnetic particles in a viscous fiuid 
l32 | , and the extensive experimental observations of 
these magnetic holes, there is a lack of theory describ- 
ing this detailed effective pair interaction potentials. No- 
tably, there has been no satisfying explanation so far for 
the existence of stable configurations of particles popu- 
lations with finite separation distances in external fields 
consisting of a circular rotating inplane component and a 
constant normal one (reported in j2]) , or for the existence 
of out-of-plane buckled structures |,8] , and no theoretical 
framework for the infiuence of the ferrofiuid layer thick- 
ness (separation of the embedding plates). 

We will show here how the magnetic boundary condi- 
tions along the confining plates lead to rich effective in- 
teraction potentials rendering for those structures, rather 
than the qualitative magneto-hydrodynamic effects pro- 
posed in j2|. In particular, we give an explanation for the 
existence of a finite equilibrium separation between par- 
ticles. Theoretical work has already been done along this 
line [Hi, but in a reduced case of constant normal field. 
The present study includes a circular high frequency os- 
cillating field in addition. The potential derived should 
be an essential brick in all the applications mentioned 
above of the magnetic holes, and in general this type of 
contribution of the confining structure should be relevant 
to any quasi-2D colloidal system with a significant dielec- 
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Figure 1: Pair of non magnetic particles in a ferrofluid layer 
and related effective dipolar moments due to the external 
magnetic field. 



trie or magnetic permeability contrast between the fluid 
medium and the confining structure, as in j^. 

In this paper, we first describe the system under study 
and review briefiy the basic modelling asumptions and 
standard theory. We next derive the instantaneous pair 
interaction potential with a precise care for the magnetic 
permeability contrast of the boundaries, and average it 
over the short oscillation periods to get the eflFective in- 
teractions. We turn then to the properties of the equi- 
librium configurations, sketch a simple dynamical theory 
and compare the theoretical and experimental results for 
the dynamics of a particle pair. Eventually, the three di- 
mensional aspect of the ferrofiuid layers is taken into ac- 
count, to evaluate both gravity-induced corrections, and 
the stability of the system under out-of-plane buckling. 



II. SYSTEM UNDER STUDY AND BASIC 
ASSUMPTIONS 

The system considered consists of two non-magnetic 
spheres inside a ferrofiuid which is homogeneous at their 
scale, whose susceptibility and magnetic permeability are 
denoted respectively Xf ^ind /z/ = -\- Xf)i where 

/io = 47r • l{)^''H.m~^. This ferrofuid is itself embed- 
ded between two glass plates considered as perfectly 
plane, parallel and nonmagnetic. The magnetic field 
anywhere between the glass plates is then decomposed 
between a uniform order zero component resulting from 
the outer imposed field, plus a perturbation due to the 
spheres. This perturbation is essentially dipolar: an iso- 
lated sphere in a field with a given uniform (constant) 
boundary value H at infinity provokes a purely dipolar 
perturbation outside it, generated by an effective dipole 
as shown in Figure 



(T = -VxH 



1 + 2X//3 



(1) 
(2) 



where x is an effective susceptibility including a demag- 
netization factor, i.e. whose precise form results when the 
boundary conditions for the magnetic field on the sur- 
face of the s phe re are properly taken into account (see 
for example |2(1 l33 |l - SI magnetic units are adopted 
throughout this whole paper, V = 7ra'^/6 and a refer re- 
spectively to the volume and diameter of a sphere. This 
result justifys the name of magnetic holes used for those 
nonmagnetic spheres, and stays valid to leading order 



when a system more complicated than an isolated sphere 
is considered: it holds as soon as the external magnetic 
field H can be considered as uniform at the scale of a 
sphere, which will be assumed and commented further in 
section IvTbI 

The average magnetic field inside the ferrofiuid H itself 
is simply related to the external uniform magnetic field 
imposed outside the glass plates H*^ through 



H = H\ 



1 



1 



XJ 



(3) 



to fulfill the Boundary Conditions along the glass- 



ferrofiuid planar interfaces - namely, H\ 



if 1 1 and 



- where the parallel and normal com- 
ponents are oriented relative to the glass plates. 

The simplest model for a particle pair can be obtained 
first by neglecting the effect of the nonmagnetic bound- 
aries: considering a couple of two such spheres of identical 
moments <t with a separation vector from center to cen- 
ter r, the total interaction energy of the system is, after 
Bleaney and Bleaney (83 |. 
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1 - 3 cos2 9 



(4) 



where r = \\r\\ and 9 = (ffo-) is the angle between the 
field and the separation vector. 

We consider now external fields composed of a circular 
in-plane component oscillating at frequency superim- 
posed to a constant normal one. At high enough fre- 
quency, the relative displacement of the particles during 
an oscillation period is negligible compared to its aver- 
age value, and the time dependence of the separation 
vector can be decoupled in a slow and a fast varying 
mode, r{t) — r{t) + Sr{f;t), with Sr — and Sr f 
. Throughout the remainder of this paper , bold sym- 
bols refer to vectors, normal ones to their norm, and 
the upper bar refers to averaged quantities over oscillat- 
ing time. Moreover, we will now use the additional con- 
straint on r that it should be in-plane, under the effect of 
the glass plates over the microspheres which center them 
at positions equally separated from the top and bottom 
boundaries. The presence of the boundaries can indeed 
be represented by equivalent image dipoles outside the 
confining plates, as estabHshed in section ITTTI which re- 
pel the dipole from the boundaries. This constraint will 
be adressed specifically in section IVIII to show that this 
centering magnetic effect is valid in most work cases when 
no lateral confinement is exerted in the system - section 
IVII Bl . The possible sink caused by gravity under the 
density contrast between the microspheres and the fiuid 
is generally neghgible - section FVII Al Thus, in the refer- 
ence frame (f ; n(g)f ; n) - where hats refer to unit vectors, 
instantaneous fields read 

ff = ij(cos(/)sina;sin0sina,cosa)"^ (5) 

with by definition 

cotana = /? = H^/H\i = Hl/{1 + Xf)Hn, (6) 
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[r;rj, comes cos 
52 \ 



7)/(l+/3^), and 



A— (2/32 _ 1 - 3cos(20 - 27)) 
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(7) 
(8) 



Consider distant enough particles to prevent contact and 
significant hydrodynamic interactions during an oscilla- 
tion period of the external field: the neglect of inertial 
terms allows to balance magnetic interaction forces and 
Stokes drag for both particles, which gives to leading or- 
der in i5r/f, 

"in-qa ■ (F + 5r) = -V[/(r) 



3Aa^ 



[{2(3^ + 3 cos(20)F + 2 sin(2(/))n ® f] (9) 



Neglecting inertial terms is easily justified since a large 
upper bound of Reynold's number can be evaluated as 
being Re = pxa^u/rj < 10~^ where x = 5r/a is the rela- 
tive amplitude of the oscillations which will be shown 
straightforwardly to be typically below 10"^, for typ- 
ical diameters a = SO/im, ferrofluid's viscosity 77 — 
9 • 10~^ Pa. s and density p = lOOOfcg.m"'^, and field 
oscillation frequency v — 100 Hz. Random thermal mo- 
tion in the ferrofluid is essentially irrelevant at these size 
scales, as will be shown in section IVIBI justifying the 
use of deterministic dynamics instead of a Brownian one. 
The above Eqs. Ij7l9|l establish that the slow motion f is 
driven by an effective potential obtained through time- 
avera ging ov er the oscillations of the field, i.e. simply 
using cos^ (/) = 1/2 at fixed r = f , 
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with U{r) = A^{2P^ - 1) 



(10) 
(11) 



while small and quick elliptic oscillations are performed, 
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[3 sin(47ri'i)r + 2 cos(47r;/t)n ® r] 
A*o(l + X/)x'i?|f 



11527r?7 



The relative magnitude of the fast oscillations 5r/f — 
Vctt^/uf^ are indeed negligible as soon v ^ Vc — 0.\Hz 
for typical ferrofiuids, Xf — 1-9; X = 0.84, 77 = 9 • 
10~^ Pa. s and fields H\\ = 14 Oe. Vc , inverse of the 
viscous relaxation time of the system, is the critical fre- 
quency introduced in above which a particle pair 
cannot anymore follow the direction of an external rotat- 
ing field due to the fiuid drag. In this paper, we study 
regimes where v > 10 Hz, for which the relative varia- 
tions of the separation vector are below 1%, and focus on 
the slow motion of the particles r(t) driven by U{r). 

In this simple picture however, this average potential 
is a simple central one whose inverse cubic range re- 
fiects the dipole-dipole nature of the interactions, and 



in the absence of any characteristic length scale, this ba- 
sic model predicts a very simple behavior for the parti- 
cle pair: depending on the ratio (3 of the normal over 
the in-plane field, either the two particles will repel each 
other without end if (3 > f3c = 1/V2, or they will attract 
each other when (3 < f3c, until the magnetic forces are 
balanced by contact forces or very short-range hydrody- 
namic forces sensitive when particles almost touch - when 
{r — a)/a gets insignificant: this theory is obviously in- 
sufficient to render for the finite equilibrium separation 
distance, sometimes at a few diameters, which is exper- 
imentally observed for a whole range of imposed fields 
A proper treatment of the boundary conditions of 
the system along the glass plates, introducing the plate 
separation as an extra length scale to the problem, will 
in the next section be shown to remedy this problem. 



III. EFFECT OF THE BOUNDARIES ON THE 
INSTANTANEOUS INTERACTIONS 

The boundaries between the ferrofiuid and the em- 
bedding glass plates are supposed to be perfectly plane. 
The two microspheres are supposed perfectly centered 
between the glass plates, and the perturbation of the 
magnetic field due to the presence of those spheres is 
modeled as a perturbation due to two identical point-like 
dipoles (T — ~xVH located at the center of the spheres. 
To fulfill the magnetic boundary conditions - i.e. the 
continuity of H\\ and B^- along the plates, a direct use 
of the image method (e.g. Weber yal) shows that this 
magnetic perturbation between the plates is equal to the 
field emitted, in an unbounded uniform medium of sus- 
ceptibility X/, by an infinite series of dipoles: the two 
original ones, at locations defined as and a;, plus an 
infinite set of images for each of them corresponding to 
the mirror symmetry across the plane boundaries of the 
sources and all of the successive images. A magnification 
factor 



Xf 



+ Mo X/ + 2 



(12) 



multipHes the amplitude of the dipolar moments at each 
symmetry operation, i.e. explicitly defined by the follow- 
ing conditions: 

For any dipole source <t at position rg, with h the 
normal separation vector between the plates, an infinite 
set of dipolar images indexed by Z € Z is defined by their 
locations and moments - see figure [21 



ri-ro = Ih 
•^'ii = '^''''^11 



(13) 
(14) 
(15) 



The total interaction energy of such a system - where 
all dipoles do not have anymore the same moment - is 
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where ^4 is the constant defined in Eq. ijHJ. The first 
term, denoted by /g'*, is the direct interaction between 
the two sources, the next one if'* corresponds to the in- 
teractions between the sources and their own images, and 
finally the term /f° corresponds to the cross-interactions 
between a source and the images of the other one. 



Figure 2: Image dipoles representing the Boundary Condi- 
tions. 
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where the sum on the b index runs over the two source 
dipoles, and the one over the c index runs over the whole 
set of sources and image dipoles - as seen in the preceding 
section, separation vectors can be considered as constant 
over the quick variations of cr, and the upper bars over r 
are implicit for the remainder. 

We will then use the following straightforward geomet- 
rical equalities resulting from Eqs. II1I6I13I15II : if the c 
index represents the Z-th image of the source b, then 
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If on contrary c represents the l-th image of one source 
as indexed in FigureE)- , and b the other source, then 



lo2 



0-b-cr,_l + (-l)'/3' 



{(Tb-rbcXcTc-rbc) 



I + /32 

{y cos (j)+l3{~iyi){y cos (j)+l3l) 



(l + /32)(y2 + ;2) 
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where 



y = h 



(20) 

(21) 
(22) 

(23) 



is the ratio of the particle in-plane separation distance to 
the plate separation. Introducing those equalities in Eq. 



IV. TIME-AVERAGED EFFECTIVE 
INTERACTIONS 

A. Derivation of the potential 

The Jf-fields consist of a perfectly circular in-plane 
component ify = (cos 0r -I- sin (f>n A r), and a normal 
component H± which is maintained constant. During 
an oscillation of the in-plane field, the only significantly 
varying quantity in Eq. itTfill is the angle (j), with once 
again cos^ cf) — 1/2. 

The term /f*, e.g. the interactions between a source 
and its own images, naturally does not depend on the 
in-plane separation vector x, and produces no net in- 
plane force. It is worthwhile to note that this would 
however not be the case for any non-plane interfaces, giv- 
ing the possibility to quench any geometrical property of 
the roughness of the interfaces in this potential. In the 
current hypothesis of purely planar plates, this interac- 
tion term is only responsible for a normal centering force, 
as will be established in section FVII Bl For the present 
purpose where the microspheres are constrained on the 
half-plane between the plates, the in-plane forces are the 
only relevant ones, and this term is simply discarded in 
the following. 

The remaining terms produce through time-average 
the interaction energy 

I7(.)^4.(J) (27) 

with A the constant defined in Eq. Q, and a dimension- 
less term 



u{y)M2p' - l)y 
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(28) 



4)'/32 



(2/2 +;2)3/2 



3y2 + 2(-l)'/2/32 

2 (y2 + ^2)5/2 



The first term reduces to the expression of the first-order 
theory of Eq. lfTT|l . which is a test of self-consistency, 
since an infinite medium would be equivalent to the ab- 
sence of permeability contrast along the plates, k — 0. 
The following ones render for the interactions between a 
source and the images of the other one. Due to the cyHn- 
drical symmetry of the problem with a purely circular 
in-plane field, this interaction is isotropic: the depen- 
dence on the separation vector x enters only through its 
norm x = hy. 
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B. Properties of the isotropic interactions 



(a) Purely attractive, /3 < /3„ 



The introduction of those images is responsible for pos- 
sible finite separation equilibrium distances for a whole 
range of /? characterizing the imposed magnetic field, as 
is illustrated in Figure El In this example, a typical sus- 
ceptibility Xf — 1-9 was considered for the ferrofluid, i.e. 
K — 0.49. The four potentials u{y) represented corre- 
spond respectively to = 0.4, 0.58, 0.8, 2.1. They were 
obtained by truncating the sums in Eq. Ij28|l to order 
I — 10, which corresponds to a relative error lower than 
lO^'^ in u for any y. 

The general typology of those potentials can be clas- 
sified in four cases as will be demonstrated in the next 
section, separated by three particular values of (3, re- 
ferred to as (3m, Pc and - which depend only on the 
susceptibility xf- 

• (a) For /? < /3m , u is a monotonically increasing 
function, and the magnetic forces are purely attrac- 
tive, thus leading any pair of spheres to contact. 

• (b) For (3jn < (3 < Pc, u presents a short-range 
attractive core, and presents both a maximum - 
an instable equilibrium point - at some distance 
yi, typically slightly below 1, and a minimum - a 
stable equibrium point - at a greater distance ys, 
usually above 1. 

Two locally stable equilibrium configurations are 
possible in principle, depending on the ratio ya — 
a/h of the particle diameter to the plate separa- 
tion, and of the initial particle separation yinit: if 
yinit > Ui, the particles should end up at the equi- 
librium separation ys, which since ys > I > ya 
corresponds to an equilibrium configuration with- 
out contact between the particles. If on contrary 
2/init < Vi, the particles should attract each other 
and end up in contact at ya- Since any separation 
y < ya is forbidden due to contact forces, this sec- 
ond case is only possible when ya < yi, i.e. when 
the ratio of plate separation over particle diame- 
ter is sufficiently big. In that case, if the thermal 
fiuctuations are large enough to let the particles go 
over the energy maximum at yi with a significant 
probability over the observation time, only one of 
the two possible equilibrium configurations will be 
thermodynamically stable, and the other one will 
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(b) Metastable, /3,„ < /3 < f3c 



(c) Stable, f3c</3<l3u 



1 2 3 4 5 6 

(d) Purely repulsive, /3u < /3 

Figure 3: Possible types of the interactions depending on f3. 
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be only metastable. The selection of stability / 
metastability between the two is determinated by 
the comparison of u{ya) and u{ys). 

• (c) For Pc < l3 < Pu, u presents only a global 
minimum, which thus corresponds to a stable equi- 
librium separation distance Us. Since generally 
Us > 1 > Ua, this equilibrium configuration corre- 
sponds to a finite separation distance hys ~ a > 0. 

• (d) For Pu < P, u is monotonically decreasing, and 
the time averaged magnetic forces are purely repul- 
sive at any separation. 



V. EQUILIBRIUM SEPARATION DISTANCE 
AS FUNCTION OF THE APPLIED FIELD 

The separation equilibrium distance yeq - which is sta- 
ble or not - correspond to the extrema of the potential, 
and can be obtained in principle by solving 



du , 
dy 



{yeq) = 0. 



(29) 



yeq corresponds to ys or yi defined in section ITV BL ac- 
cording to the sign of d^u/dy"^ . Derivating Eq. II28|I with 
respect to the scaled separation y, leads straightforwardly 
to 



y 



^du 



+ 00 



3 dy 

My) 



(1 - 2/32)-2^«;'AKy)(-l + 2(-l)'/32) (30) 



{y' - ^l^)y' 

(y2 ^2)7/2 



(31) 



The Eq. l(3(Hl above is the sum of a term independent of /?, 
plus another proportional to The constant term can 
be shown to be positive, and the prefactor of 0^ strictly 
negative, for any possible (j/, k) - i.e. any ?; > 0, < k < 
1. Thus, du/dy{y, (3) is a monotonic decreasing function 
of , equal to zero when 



Poiy) = 



i + 2Er=i^'My) 

2 + 4EtT(-l)'«^'AK2/) 



(32) 



Thus, for a given field configuration and separation 
y, pair interactions are attractive, i.e. du/dy > 0, if 
f3 < Poiy), and conversely if /? > f3o{y). A numerical 
study of the above function, for any k, shows that Po{y) 
is monotonically decreasing from /3c = l/\/2 to a finite 
positive minimum /?„('«), between y = to ym^n), and 
next monotonically increasing up to a finite limit /?«(«;) > 
f3c between ?/m(K) and y -\-oo. 

These considerations allow to obtain by a direct graph- 
ical inversion of /3o(y), the possible roots ysiP) and yi{l3) 
for which the interaction forces are zero at a given field 
geometry (3, as shown in Figure which is obtained for 
the particular case Xf = 1-9: i-S- k = 0.49. This nonlin- 
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Figure 4: Equilibrium separation of a pair as function of the 
applied field 



ear dependence of the equilibrium separation ysiP) seems 
more complex than observed in earlier experiments by 
Helgesen and Skjeltorp Q. This apparent discrepancy 
will be resolved in the next section, which is centered on 
finite time results. 

Identifying the three parameters (}m,Pc,Pu identified 
above as 



Pc = Po{y = o) 

P,ri = min/3o(y) = /?o(2/m) 
y 

Pu = lim Po{y), 



(33) 
(34) 

(35) 



these arguments prove that the pair effective potentials 
belong to one of the four types described in the previous 
section: 

(a) \i P < Pm, the potential is purely attractive at any 
separation. 

(b) \l Pm < P < Pc, there are two roots to Eq. ll29)l . 
denoted yi{P) and ys{P)- the potential is attractive below 
yi or above ys, and repulsive between both. 

(c) If Pc < P < Pu, the potential presents a single 
minimum at ysiP)- 

(d) If P > Pu- the potential is purely repulsive and 
there is no equilibrium separation. 

The definition of Pu above - Eq. (j35|l - shows also 
clearly that lim^^^-(j/s) = -t-00: in principle, it should 
be possible to drive a pair of microspheres in an equi- 
librium configuration with any desirable separation dis- 
tance. Naturally, since the magnetic interactions decay 
rapidly with distance, thermal processes or any kind of 
external perturbation in the fiuid fiow, or default in the 
planarity of the plates, will be predominent at large sep- 
arations, where this theory will become inapplicable. 

The dependence of Pc,Pu,Pm on the susceptibility of 
the ferrofiuid (through the parameter k) is as follows: 
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and similarly for 
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Figure 5: Equilibrium separation distances as function of the 
applied field, for various ferrofluids. 



Replacing = in Eq. l(32)l shows that 

Pc = I/V2 (36) 

independently of k. Similarly, since liniy^+oo ^i{y) — 1, 
/?„ is easily summed as 



/3u 



V2Vi+2ES(-'^y 

1 l + K 



(37) 



A numerical study of /3m shows that it decreases mono- 
tonically with k, down to zero when k ^ 1. 

This shows that the range of (3 over which stable equi- 
librium distances exist is larger when the susceptibility 
of the ferrofiuid is important (k increases with Xf^Pm 
and Pu are respectively decreasing and increasing with 
k) , up to the limiting case of an infinitly susceptible fer- 
rofiuid K ^ 1 {xf » 1), for which f3m ^ 0, /?„ ^ +00, 
and there is a finite stable equilibrium separation for any 
ratio p. 

The other limiting case k ^ is obtained directly by 
discarding the sums in Eq. Ij32|l - their convergence to 
in that limit for any y is straightforward -. This shows 
that limK^o Pm = limK^o /?« = /3c = 1/ V^'- without any 
permeability contrast along the glass plates, no images 
are felt and the first-order theory is recovered - the in- 
teractions are simply purely attractive when /3 < /3c, and 
purely repulsive when (3 > f3c, with the stable regimes 
/3 e [/3m, /3u] disappearing. 

An overall picture of yeqiP) for different values of k 
distributed regularly between and 0.9 is given in figure 
©. The sums in Eq. II32II have been truncated to or- 
der 10, which results in an accuracy better than 1% for 
the displayed function /3o (y) - functions A/ (y) defined in 
Eq. (Ell can be bounded by \Xi{y)\ < l/(y^ + /^)^/^ so 



that 



Thus, neglecting 



terms from order /V = 11 in the sums, for any k < 
0.9, results in a relative error on (3o{y), smaller than 
0.9"/(l - 0.9)11^ ~ 0.0023. 



VI. FINITE-TIME THEORY AND 
EXPERIMENTS 

A. Simple time-dependent theory 

The preceding section was centered on equilibrium 
properties of this system, but the relaxation time to reach 
equilibrium amounts to hours or days in certain configu- 
rations, as will be shown here. In order to compare effi- 
ciently theoretical and experimental data, we have there- 
fore concentrated on the slow dynamics of this system, 
starting from hole pairs in contact under the effect of a 
purely in-plane field /3 = 0, through the following scheme: 
neglecting once again the inertial terms. Stokes drag and 
magnetic interactions are balanced to obtain 



dx d — 
STTna—— = — —Uix) 
' dt dx ^ ' 



(38) 



where the effective potential includes the images due to 
the boundaries, Eq. Ij27|l . The viscosity 77 above is renor- 
malized to take into account hydrodynamic interactions 
with the confining plates, as will be detailed further. 
In dimensionless units, this equation leads to 



T 



dy_ 
dt' 

with T 



dy 



8647] 



/i5 



Mo(l + X/)x'i/|f 



(39) 
(40) 



and y{t') = x/h, t' — t/T, and u{y) the potential defined 
in Eq. Ij28|l . The characteristic time above lies typically 
around 30 s to 5 min at usual working parameters, as will 
be shown in the next section, and moreover the potential 
wells can be pretty fiat, thus producing often metastable 
situations that last from minutes to hours (the driving 
force close to equilibrium position is proportional to the 
distance to it, times the second derivative of the potential 
in the well, and u"{ys) when (3 ^ j3~). 

Starting from a pair configuration in contact, and set- 
ting at time the field parameters (/3 ratio and mag- 
nitude i/||) to a constant value, the time t'{y) to reach 
a given separation will be directly obtained through nu- 
merical integration of the differential Eq. (j39ll : 



t'iv) = -T 



-a/h 



dz/u{z). 



(41) 



An inverse representation of the above is plotted for a 
ferrofiuid of susceptibility Xf = 1-9- In FigureEl the solid 
thin line represents y{t') at fixed /3 = 0.8, and in Figure 
the dashed lines represent y{f3) for four characteristic 
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times t' = .25, 1, 4, 16, slowly converging to equilibrium 
separation ys{P), plotted as a solid line. 

To compare this with experiments, the fine tuning of 
the time dependence requires a refined analysis of the hy- 
drodynamic interactions in the above: since experiments 
are carried out in cells of width h comparable with the 
diameter a of the embedded holes - typically h/a lies 
between 1.1 and 2 -, a strong hydrodynamic coupling 
with the confining plates is present. Considering that 
both particles sit at a fixed fraction z of the plate sep- 
aration relative to the central position (nonzero z can 
be obtained in principle between horizontal plates due to 
the density contrast between the particles and ferrofiuid) , 
these interactions are represented according to , by a 
normalization of the Stokes drag as 

v/m = f[a/ih + 2zh)] + f[a/ih - 2zh)] - 1 (42) 

with r]o the naked viscosity of the carrier ferrofiuid, and 



fix) =1 \ 

' ^ 16 8 256 16 



+ 0{x^). (43) 



Two experimental cases will be considered here, where 
h/a = 1.4; z ~ and h/a = 2, 2 ~ 0.14, corresponding 
respectively to rj/riQ ~ 2.6 and 2.4. This is consistent 
with Faucheux and Libchaber's measures for a confined 
brownian motion , and with an experimental value 
77/770 — 2.4, that we measured in the h/a — 1.4 case by 
placing 50 iim diameter particles between 70 distant 
plates set up in vertical position, without any magnetic 
field, and by recording the motion of a single particle 
under the effect of buoyancy forces. The density of the 
ferrofiuid was pf — 1.24, the one of the polystyrene par- 
ticles pp ~ \. This value is slightly below the theoret- 
ical 77/770 ~ 2.6, which is consistent with the effect of 
the brownian motion of the particle along the normal 
direction, as analyzed in We have chosen to use 

77/770 = 2.5 in the following. 

Hydrodynamic interactions between both particles 
were neglected, which should be a relatively poor ap- 
proximation for particles close to contact, but become 
reasonable at separations x/h > 2 where most of the 
time is spent to achieve equilibrium at relatively large 
separations, and is therefore the most important one in 
the present context. The typical magnitude of this error 
can be roughly estimated through the analyzes performed 
by Dufresne et al. or Grier and Behrens (21I |^ for a 
close context: they derived for two particles of diameter 
a, separated by x and at a distance h/2 from a single 
plate, the hydrodynamic corrections to the mobility to 
first order in x/a and x/h. Considering for simplicity a 
double contribution for two plates relatively to Eq. (13) 
in jiJi , contributions due to the particle-particle interac- 
tions, become for h/a = 1.4; x/h > 2, less than 30% of 
the one due to the plates - term 9a;/16 in Eq. Il43|l . 

We neglect also the rotational degrees of freedom of 
the ferrofiuid itself, which can lead to the rotation of the 



nonmagnetic spheres, and induce another type hydrody- 
namic interactions between pairs of spheres: in ferrofiu- 
ids submitted to circular magnetic fields, the rotation of 
the magnetite particles induces asymmetric stresses in 
the fiuid, which leads to a counter-rotation of the mag- 
netic holes 0,123 ■ The mismatch between this rotational 
motion of the magnetic holes and the one of the mag- 
netites could in principle induce a vortex in the fiuid 
fiow around each of the holes, which would lead to a net 
hydrodynamic torque over close pairs of holes. Nonethe- 
less, with the ferrofiuid and field frequency regime used 
here (<100 Hz), this rotational motion of the holes is 
so slow that it is hardly detectable with the field fre- 
quency used {i^f < lOOHz). Experimentally, with a sim- 
ilar ferrofiuid (kerosene-based, x ~ 0-8), the hole's fre- 
quency was bounded by < O.OlHz. This rotational 
motion was described theoretically in detail by Miguel 
and Rubi Through this theory, for the ferrofiuid 

and field frequencies used here, the frequency of the mag- 
netic holes will straightforwardly be shown to be lower 
than Vs < 0.05Hz. This justifies for the present study 
the neglect of these rotation-induced hydrodynamic in- 
teractions. These interactions would in any case lead 
to a purely rotational motion of the hole pair, decoupled 
from the purely central forces induced by the magnetic in- 
teractions studied here. Experimentally, some very slow 
rotational motion of the hole pairs was indeed occasion- 
ally observed when the particles were close (non periodic, 
angular velocity always lower than 0.001 Hz), but no sys- 
tematic trend for its direction or velocity was noted, i.e. 
this effect was beyond the experimental error. 

To obtain the theoretical estimate of the hole's fre- 
quency above, the Langevin parameter of the ferrofiuid, 
defined f3| as = rnQH/ksT where ttio is the magnetic 
moment of the magnetites, is derived from the satura- 
tion magnetization of the ferrofiuid and its susceptibility 
[3fll | as p, — 3H/xMsnt. For the ferrofiuid used, x — 1-9, 
Msat = 200G, = 14 Oe, and p = 0.26. This enables to 
express the ratio of the hole's frequency over the field's 
one, as 

Vs/vf = -^{p - tanh^)/(^ + tanh^) = -0.00041 (44) 

where $ = —0.036 is the magnetite volume fraction of 
the ferrofiuid. For the regime vf < lOOHz where the 
experiments are carried, the spheres frequency is thus 
below i^s = 0.041iJz. 



B. Experimental results and scaling 

We used pairs of a = 50 pm diameter neutral 
polystyrene spheres of density Pp = 1, designed accord- 
ing to Ugelstad's technique ^4Qj, and ferrofiuids of sus- 
ceptibility Xf = 1.9, density pf = 1.24 and viscos- 
ity Vn — 0.009 Pa. s [4l|, confined in horizontal cells of 
thickness 70 pm or 100 pm and width of order centime- 
ters. The thickness of the cell was obtained by confin- 
ing plates by quenching a few 70 pm diameter spheres 
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between the two plates clamped together (these spac- 
ers where typically half a centimeter distant from each 
other) . The oscillating inplane and constant normal field 
were generated by external coils, with a typical magni- 
tude = Ho — 14.2 Oe, with frequencies from 10 to 
100 Hz. Magnitudes and phase of the field are accurate 
up to 1%, which also corresponds to the degree of ho- 
mogeneity of the in and out-of-plane fields throughout 
the entire cell. The direction of the constant field varies 
slightly along the cell, with a maximum 3° misalignment 
from the direction normal to the confining planes (these 
accuracies for the homogeneity and misalignment of the 
field were obtained directly by considering the geometry 
of the Helmholtz coils generating the field, whose charac- 
teristic extent is 10 cm, together with a .5 mm accuracy 
for the position of the cell inside them). The in-plane 
motion of the particles was recorded using a microscope, 
and a camera linked to a numerical data analysis setup. 
The entire experimental setup is described, for example, 
in Ref. The samples were prepared with a pure 

inplane field to bring the particles in contact, after that 
the inplane field was maintained constant and the nor- 
mal field was set to a constant H^ = (1 -I- x/)/3i/||. The 
normal field required typically a few seconds to stabilize. 
Separation between the particles were then recorded ev- 
ery 10 s, for 30 mm. The concentration of the spheres 
in the entire sample was such as the nearest sphere or 
spacer would sit at least 20 diameters apart from the ob- 
served pair, which was sufficiently dilute not to infiunce 
the motion of the observed pair. 



Time-dependent result at fixed field geometry and 
discussion. 



The scaled separation as function of time is shown in 
Figure El for five experiments carried out at /? = 0.8, 
i.e. for the potential represented in Figure 13 (c) , which 
presents a single minimum at ys(0.8) = 2.35. Four ex- 
periments were carried out in a /i = 70 /im thick cell, two 
of them at identical field amplitudes iJy = Hq — 14.2 Oe, 
and two other at = 0.7 Hq; 0.5 Hq. The last experi- 
ment was carried out in a /i = 100 thick cell, with an 
inplane field Hq. These parameters corresponded respec- 
tively to characteristic times evaluated using Eq. Ij4()|l as 
T = 32 s, 63 s, 129s, 192s. 

a. Comparison with the theory. The large part of 
the figure represents the short time evolution t <T (typ- 
ically the first minutes), the encapsulated part shows 
longer time (typically 30 min). Every data point (sep- 
arated by 10 s) was used for the early regime, one point 
out of 2 to 12 depending on the experiment for the later 
one. At shorter times, the scaled data collapse reason- 
ably well on the theoretical curve obtained from the the- 
ory sketched so far - solid line. The only two main out- 
liers (first filled circle and triangle, t = 10s) presumably 
corresponded to a weaker normal field in the first few sec- 
onds, until it stabilized. For illustration, the first-order 
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Figure 6: Scaled separation as function of time for various 
field magnitudes and plate separations, for /3 = H^/H\\ — 0.8. 
The inset represents the same data for longer times. 



theory neglecting the magnetic confinement leads to the 
dot-dashed curve, obtained using Eq. Ij41|l with a bare 
potential retaining only the source-source term. This too 
simple theory clearly differs from data both at short and 
large times. 

h. Corrections due to buoyancy. A third dashed the- 
oretical curve was plotted, corresponding to a slightly 
refined calculation of the interaction potential, where 
buoyancy forces due to the density contrast between 
ferrofiuid and particles led to a shift z of the particle 
pairs from the central mid-plane of the cell. The rela- 
tive displacement z/h due to this gravitational correc- 
tion is evaluated in section IVII Al as 15% in the weakest 
fields 0.5 iJo, h/a = 1.4 case, and in the thicker layer 
Hq, h/a — 2 case, whereas it should remain around 6% 
for the thiner cells, higher field Hq case. The extended 
potential to take this lateral shift into account is also de- 
rived in section FVII AL and the dashed curve corresponds 
to this potential at a fixed shift z/h = 14%, which is 
the maximum possible in the thin cells, since it would 
correspond to contact of the particles with one of the 
plates, coinciding with z/h = (1 — a/h)/2 at h/a = 1.4. 
Considering this correction, one would expect the dia- 
mond and squares to follow the dashed line, the circles 
to be close to the solid line, and the triangle to sit in be- 
tween. This is indeed approximately the case at shorter 
times for the diamonds and rectangles, and the filled cir- 
cles fit well with the centered theory, both at short and 
long times. Nonetheless, it is worth noticing that the 
open circles, collapse rather with the lower field exper- 
iments than with the filled circles, another experiment 
carried with identical parameters. This gives an estimate 
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of the reproducibility of these experiments, correspond- 
ing roughly to a relative experimental error bar of 10% 
for the separation at a given time, estimated between 
open and filled circle cases. The gravity- induced correc- 
tion discussed above is of order 5% for the separation as 
function of the dimensionless time, and the effect of this 
shift on the magnetic interactions can therefore hardly 
be distinghished from experimental dispersion in terms 
of separation. Still, the renormalization of the time due 
to hydrodynamic interactions with the plates would lead 
in the thick cell case h/a = 2, to Ty/r/o — 1-2 if the parti- 
cles were considered as centered z/h = 0, instead of the 
value ri/rjo = 2.5 that we used and corresponded to the 
predicted shift in that case, h/a — 2, z/h — 0.15. This 
gravity-induced correction is then clearly sensible for the 
hydrodynamic corrections, if not so much on the mag- 
netic interactions, since the above incorrect viscosity i] 
would correspond to dividing T by 2, which would dou- 
ble the abscissas of the square data points and drive them 
way out of the theory and rest of the data. This shows 
qualitatively that this shift was indeed present in cases 
h/a = 2, Hq, and h/a = 1.4, 0.5 Hq, and that particles 
sat close to or in contact with a plate in this last case. 

c. Corrections due to experimental error on the field 
directions. The fourth dotted curve represents the theo- 
retical effect of a misalignment of angle a = 2.5° between 
the constant field and the normal direction. A direct gen- 
eralization of sections^tol^to such configurations with 
a slightly tilted constant field, shows that the time aver- 
aged potential is still of the form in Eqs. II25I26|I . with a 
modified parameter 

P' = ^1 

(1 + x/)(J?|| + V2am{a) cos{0)Hl) 

~ /3- V2(l + x/)/3^sin(a)cos(6l), (45) 

where is the angle between the projection of the con- 
stant field over the plane and the separation vector. This 
modified interaction potential leads to a torque tending 
to align the particle pair with the direction 9 = 0, and a 
radial interaction force corresponding to a modified P' in 
[(3 — V2{1 + Xf)0^ sin(a); /?], (}' decreasing with time to- 
wards the lower limit as the pair aligns with the inplane 
constant component of the field. The dotted curve cor- 
responds to this lower limit for a possible misalignment 
a = 2.5°, which is evaluated from the above as /?' = 0.68. 

For the longer time period shown in the encapsulated 
part of figure EI an apparent equilibrium position was 
reached in each experiment after typically 6 T - no more 
than a 1% relative motion was noticed later when the 
experiments were several hours. This equilibrium sepa- 
ration corresponds theoretically to the equilibrium one 
studied in the previous section, ?/s(0.8) = 2.35. The 
high field and thin cell experiment (filled circles) agree 
well (within 2%) with the theory with a purely normal 
field (soHd line), but discrepancies between soHd line and 
experiments are noticeable at longer times for the four 
other experiments. A comparison of the data with the 



dotted line shows that a misalignment of order 2.5° be- 
tween constant field and normal direction is sufficient to 
explain these discrepancies: in these experiments, the 
particle pairs started at a relatively large angle 9 from 
the inplane component of the constant field, which is 
why the unmodified theory and experimental data are 
close for short times. At longer times, the particle pairs 
aligned with the direction = 0, and the modified theory 
P' = 0.68 agrees well with the data. An initial rotation 
of the particle pair and subsequent locking of this di- 
rection in a particular one was indeed observed in these 
experiments. 

Brownian motion in the ferrofluid can be proved to 
be entirely negligible for the relatively large particles 
and field we worked with, its relative magnitude com- 
pared to magnetic interaction energy being kT/ min U < 
kT/[A{d^/h^)mmu\ ~ lO-^/minw < 5 • 10"^ for the 
fields, particles and plate separations considered here. 



S. Scaled time-dependent results at various field geometries 

Experiments were carried out with inplane field magni- 
tudes Hq, and normal fields jumping at initial time from 
to /3(1 -I- Xf)Ho with various P from to 3.5. This was 
done for a = 50 /im diameter particles, and plates sepa- 
ration h = 70 and 100 ^m. The results, scaled sepa- 
ration as function of P, at four values of the scaled time, 
are shown in Figure E{a) for the thinner cell. The error 
bars correspond to a possible misalignment a = 2.5° be- 
tween the constant field and the normal direction: they 
represent the Hmits [P — V2{l+Xf)P'^ sin(a); /?] for the ef- 
fective P' parameter as explained in the previous section. 
For the "c"-regime P > Pc = l/\/2 ^ 0.7, theory and ex- 
periment agree well for any of the tested field parameters 
and times. 

The solid line represents the theoretical equihbrium 
value, studied in section ^ This is reached within typi- 
cally 16 T (8 min for h/a = 1.4, = Hq) for ^ < 1, or 
longer time at higher p. This is the main reason why the 
upwards curvature of the theoretical solid curve at larger 
values of P is not observed in experimental data, which 
correspond to finite times, and for which other types of 
perturbations always enter the picture at very large times 
and distances. 

In Figure Cfb) , we present the results of experiments 
carried out at two different plate separations, as func- 
tion of time and value of p. The error bars have been 
omitted for readability, and the experimental points rep- 
resented correspond to a constant field supposed purely 
normal (i.e. the abcissas are the upper limit of the er- 
ror bars in Figure [3(a)). The experiments carried out 
in thicker cells, corresponding to weaker magnetic inter- 
action forces, are more sensitive any perturbations. The 
relative data collapse for both plate separations at 0.25 T 
and T when (3 > Pc nonetheless show that the separations 
in this regime scale with plate separation, and not with 
particle diameter. Apart from the misalignment of the 
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Figure 7: Separation as function of time and /? for various 
field configurations, times and cell sizes: merged data and 
theory. 



constant field with the normal direction, a possible source 
for these perturbations is as follows: when the particles 
come close to equilibrium, inplane magnetic forces tend 
to zero, and the particle motion becomes more sensitive 
to any interactions with the local environment (confin- 
ing plates) - this being of course even more the case for 
weaker fields or larger hja: there seems to be a pinning 
(friction) of the particles to an absolute plate position at 
large times. The physical origin of this pinning is possibly 



due to roughness of the plates (especially when particles 
are almost in contact), which can quench the particles 
through the magnetic perturbation due to this rough- 
ness (the repulsion effect of a dipole by its images would 
make a particle sit preferably in positions of larger plate 
separation), or alternatively when particles are almost 
in contact with the plates can result in an inplane com- 
ponent of the hydrodynamic coupling or contact forces 
responding to buoyancy forces. Instead of plate rough- 
ness, the same type of qualitative effects could be due 
to small impurities in the ferrofiuids, starting to stick to 
the plates or particles at large times, when the chemi- 
cal surfactant layers around large particles and possible 
impurities, break apart in some points. 

For the regime /3 < /3c, the simple theory presented 
here would predict that particles stay always in contact 
sX y = a/d = 0.5 or 0.71 for (3 < (3m ~ 0.55, and for 
Pm < (3 < Pc would either stay in contact if a/ft. < yi{(3) 
or go to the secondary minimum ysiP) if yiiP) < a/h 
(the separation between first and second case happening 
at^ = 0.61 for a/ft ^ 1/1.4, and /3 = 0.67 for a/ft = 1/2). 
Particles seem indeed to be in contact for (3 < 0.2, but 
start to separate significantly well before Pm- We note 
also that this separation seems grossly to be proportional 
to the particle diameter when /3 < 0.3, where some finite 
separation can already be observed - ordinates of opened 
and filled symbols are multiple of each other through a 
factor 100/70 - and proportional to plate separation in 
the regime Pm < (3 < (3c- This shows that an extra 
physical effect that was not taken into account here, gen- 
erated repulsive forces, whose range is finite but scales 
with the particle diameter. This effect suppresses then 
the short range attraction in the regime (3m < P < (3c, 
so that the particle jump directly to ys{P), which is al- 
ways a stable minimum and not a metastable one. For 
P < Pm, this extra effect starts to separate the particles 
in proportion to particle diameter. The physical origin 
of this short range repulsion should not be the particle- 
particle hydrodynamic interactions, which should slow 
down the relative motion rather than result in a net re- 
pulsion - see [13, E3| A probable candidate for this 
repulsion is rather the magnetic effect of the finite size of 
the spheres, particularly sensible when particles are close 
to contact: even if an isolated nonmagnetic sphere gen- 
erates a purely dipolar perturbation when it is isolated 
in a homogeneous susceptible medium, this dipolar per- 
turbation does not fulfill the boundary conditions along 
the surface of another magnetic sphere, sufficiently close 
of the first one to feel the heterogeneity of the perturba- 
tion at the scale of its diameter, which is naturally the 
case at a finite separation/diameter ratio. To model this 
short range repulsion requires accounting for the mag- 
netic perturbation generated by this non pointlike char- 
acter of close enough spherical particles. Though this can 
be performed by a simple image method for pairs of discs 
in 2D, one can show that this does not extend to pairs of 
spheres in 3D, and the proper mathematical description 
of this perturbation requires the use of a series of spher- 
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ical harmonics, which was not performed in the present 
study. To conclude this discussion, we note that this cor- 
rection seems to be neghgible at separations exceeding 
the sphere diameter xja > 2, as shows the agreement 
between experimental data and the present theory when 
/3 > A„. 

Finally, we note that the present results are not contra- 
dictory with similar experiments carried out in |3| with a 
sHghtly different ferrofluid, plate separation and particle 
size, where it was reported that the equilibrium parti- 
cle separation is approximately linear in /3 once particles 
have started to separate: for example, this is also the 
case with the present ferrofluid, in the particular case 
hja = 1.4 in the regime 0.3 < /? < 1 - see the filled 
triangles or theoretical curve in Figure H^a). This Hnear 
property is however shown here to be a mere coincidence, 
for this does not hold in the same (3 regime for hja = 2, 
or for any hja when /9 > 1, where ysifi) is curved up- 
wards, and any result y{f3,t) observed at a given finite 
time t is curved downwards. 



VII. GENERALIZATION TO QUASI2D 
SYSTEMS 



A. Gravity induced corrections 

Although the effect of the dipolar images in the con- 
fining plates tends to center particles at a mid-plane po- 
sition, the density contrast between the ferrofluid and 
particles tends to drive the particle out of it for large 
enough particles. An estimation of this effect can be ob- 
tained by considering for each particle, the sole effect of 
its own images, plus buoyancy forces - for the simple es- 
timation we look at here, we will neglect the coupling 
between one source and the images of the other. Extend- 
ing the analysis performed in sections ITTTI and ITV Al to a 
single dipole lying at a vertical distance z from the center 
between two horizontal plates, comes directly 



U (z) = A — U' 
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This potential can be shown to be always centering for 
any value of /?, i.e. to have a single minimum in s = 0, 
and to diverge to infinity at s = ±0.5 - plate contact for 
very small particles. Equilibrium between gravity forces 
and magnetic interactions between the dipole and its im- 
ages leads to 
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Figure 8: Shifted hole pairs considered. 



For small separations (i.e, small particles or strong 
enough fields), a Taylor expansion to first order around 
the plates' center gives 



^ = 96 • C(k) • (1 + 2/32) . |s| 
ds 



E 



At2"+V(2n+1)^ 



(50) 
(51) 



which is valid up to 25% for s < 0.1. Thus, the displace- 
ment looked for can be evaluated as 



s = z/a = 



(Pg - Pf)9h 



2C(k)(1 



(52) 



when this quantity does not exceed 0.1, or directly using 
Eqs. (|49l5n|l otherwise. For the ferrofiuid and parti- 
cles we used, this led respectively for h = 7Q fim, H\\ — 
Ho;0.7Ho or 0.5Ho, and h = lOO^m; = Hq, to 
s = 0.06, 0.10, 0.15 and 0.15 . 

The magnetic interaction potentials were unaffected up 
to 1% in the s = 0.06 vertical shift case, and the dashed 
curve in Figure|2|was obtained by considering particles at 
a fixed s = 0.14 out-of-midplane shift, using a generalized 
potential obtained for a configuration sketched in Figure 
IHl through an extension of the method used in sections 
ITTTI and ITVAT as 
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(54) 



This s = 0.14 value was picked to represent the magnetic 
effect of a shift sufficient to bring particles in contact with 
the plates in the hja = 1.4; 0.5 Hq case. To an accuracy 
of 1%, the results for s = 0.15 were very close to this 
case, the ones for s — 0.06 very close to pure inplane 
situations, and the situation s = 0.10 fell roughly halfway 
between both, and were therefore omitted from Figure 
for readability. 



B. Stability of the plane solutions and buckled 
configurations 

When particles are sufficiently small or fields suffi- 
ciently high, the preceding section estabHshes that the 
confining plates have an effective repulsive effect on an 
isolated particle, which is therefore centered on midplane. 
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Figure 9: Tilted hole pairs considered. 



In the case of a pair of particles, the interactions be- 
tween one source and the images of the other one might 
nonetheless modify that picture, and make the plane so- 
lutions described in this paper unstable as have been 
observed in some experiments. Neglecting gravity, we 
will here generalize the interaction potential, to configu- 
rations where particle pairs are allowed to tilt on both 
sides of the midplane, i.e. where both are displaced 
by the same distance z on both sides of it - cf Figure 
El We consider only symmetric situations, due to the 
symmetry of the problem under parity in the absence 
of gravitational forces. The effective interaction poten- 
tial is obtained similarly to the plane one, and comes as 
C/(x, z) — A{a^ /h^) ■ u{x/h,z/h) where 



+00 

i(2/,s) = 2 5] J'l 
— 00 



1 + (-l)'/3' 



(y2 + (/+p,s))2)3/2 
2 (y2 + (/+p,s)2)5/2 



Pi 



-f2^Kl'l(l-2(-l)'/?2 

lei' 
1 + (-!)' 



(55) 



1 



\l + {2-pi)s\^ \lt 
(56) 



which reduces to the previous in-plane solution, Eq. ll2? 
when s = 0. Contour plots of this pair potential are dis- 
played in FigureEH for the ferrofluid used here, Xf = 1-9) 
and the four types of potentials - (3 values are identical 
to the ones adopted in section IIVBI Black and white 
represent respectively u = —0.2 and 0.5 in subfigures (a) 
and (b), tt = —0.2 and 1 in (c) and (d), and the grey level 
linear in u in between. Out-of-plane values represented 
cover the whole possible range < s < 0.5 in (c) and 
(d), and are restricted to 10% from the midplane in (a) 
and (b). The inplane configurations correspond to the 
bottom axis of those graphs. 

Inplane solutions are in principle locally stable if 
d^u/ds^{y,0) > 0, otherwise particle pairs will tend 
to tilt. For both first cases, /S < /?c, we note that 
d^u/ds"^ > for any possible (j/, s), and any tilt is re- 
stored by the magnetic image effect: the plane configu- 
rations are indeed stable. As soon as f3 > f3c, besides the 
minima (y, s) — (y^, 0) or (-f 00, 0) in c or d case, another 
local minimum appears at a certain (0,Se(/3)): particles 
can be stable at a finite distance on top of each other - 
the pair tends to align with the large normal field -, or 
if they are close enough to be attracted by this potential 
minimum but too large (a > Se(/5) and a > h/2), con- 
tact forces between them and with the confining plates 
will attract them to a "buckled" configuration where both 
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(a) /3 < /3„ 




O.b 1 1.5 2 2.5 3 



(b) |3rr^<f3<|3c 




(c) Pc<l3</3u 




(d) Pu<f3 



Figure 10: Contour plots of the interaction potential as func- 
tion of inplane separation y and tilt coordinate s. 
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Figure 11: Stability limits of the in-plane configurations. 

particles are in contact with each other, and with one dif- 
ferent plate each. Note that < Se(/3) < 0.5, i.e. very 
small particles a ^ at this second minimum would sit 
on top of each other, neither in contact between them nor 
with the plates. The criterion to determine whether par- 
ticles are attracted by the inplane solution, or the buck- 
led configurations, is to determine whether the present 
{y, s) lie in the basin of attraction of one minimum or 
the other. Both basins of attraction are separated by 
a ridge of the potential, on which u decreases under the 
effect of any perturbation of (y, s) apart from the ones di- 
rected exactly along its gradient. This boundary between 
both basins of attraction, noted yi{s,(3) was determined 
numerically and plotted as the grey solid line in Figure 
El (c) and (d). We note that the whole axis (y, 0) lies 
in the basin of attraction of the plane minimum, so any 
particle pair starting with no tilt should end up so. How- 
ever, at short enough distances y in a c-field, or at any 
inplane distance in a d-field, certain configurations with 
a finite tilt are attracted by the normal-aligned pair min- 
imum, and will end up in a buckled contact configuration 
or normal-aligned pair. We have determined for every (3, 
the maximum ym{P) — maxj yi{s, (3): for a given /3, when 
y > VmiP) any tilted configuration will be attracted by 
the inplane solution, whereas for certain finite tilts at 
close enough y < ym{P), the pairs will be attracted to- 
wards buckled in-contact configurations. The function 
ym{fi) is displayed as the dashed curve in Figure ITTl - 
the solid curve in a reminder of the equilibrium inplane 
distances ys{/3),yi{l3) determined in sectional The func- 
tion ymi(3) diverges when /? ^ /3~, illustrating the fact 
that in any d-type field, configurations with out-of-plane 
tilts s close to 0.5, i.e. both particle centers almost along 
the plates (for small enough particles) will be attracted 
by the normal-aligned pair mode. 

This effect is believed to be responsible for the lattices 
of buckled chains of particles in contact observed in Q. 
These were observed in pure normal fields (/? = -l-oo), 
and the non-trivial character of the lattices, being hexag- 
onal or square instead of triangular lattices characteristic 



of attractive interactions at any range, can be qualita- 
tively explained by frustration effects: neighboring par- 
ticles tend to sit in top-down contact configurations, but 
top-top or down-down configurations are repulsive - cf 
shifted potential developed in the previous section - and 
therefore lattices with noneven number of particles along 
the loops, as the triangular one, are disfavored in com- 
parison with those involving even numbers in loops, as 
the square or hexagonal ones. 

Eventually, the local stability of a plane configuration 
was investigated in fields of c- or d-type: the second 
derivative d'^u/ds^ is positive at s = for large enough 
distances, and small out-of-plane tilts will be restored by 
magnetic interactions, but below a certain y < yp{P), 
d^u/ds'^{y, s = 0) < 0, and the inplane solution will be 
locally instable. However, the preceding estabHshed that 
this out-of-plane character will be transient, since this 
case will still be attracted at long times by the inplane 
solution. This upper limit yp{f3) below which inplane 
configurations will go to a transient out-of-plane regime, 
was determined numerically for any /3 and plotted as the 
dash-dot curve in Figure ITTl 



VIII. CONCLUSIONS 

For pairs of magnetic holes in ferrofiuid layers of fi- 
nite thickness, exposed to fast oscillating conic magnetic 
fields, we have established the effective interaction po- 
tential driving their slow motion. The importance of the 
magnetic permeability contrast between ferrofiuid and 
confining plates on these interaction was demonstrated, 
and the resulting pair potential analytically derived. This 
allowed the classification of those interaction potentials 
into four types, two of which present a secondary min- 
imum at a finite distance. A simple finite time theory 
for these non-Brownian microspheres, including the hy- 
drodynamic interactions of the holes with the confining 
plates, was directly compared and confirmed by experi- 
mental results, through data collapse of the scaled sepa- 
ration as function of scaled time, for various plate separa- 
tions and field magnitudes. The relaxation time of these 
systems to reach equilibrium when there is any, was typ- 
ically a few minutes or above. Eventually, we generalized 
the study to full three dimensions in the layer thickness, 
and studied the stability of the in-plane configurations. 
This established that although inplane configurations are 
stable at small normal fields, or at large ones and suffi- 
cient particle separation, hole pairs can be attracted by 
another stable configuration of tilted pairs of particles 
with contact between them and contact with one plate 
each. 

This simple theory renders for so far unexplained ob- 
served configurations of magnetic holes, namely 2D lat- 
tices with finite separation, or buckled lattices of tilted 
pairs of particles. In principle, the effect described here 
in detail should be important in any colloidal system con- 
fined in a layer with significant magnetic permeability or 
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dielectric contrast between the fluid and conflning struc- 
ture. 

The abihty to tune the equilibrium distance at will 
through the ratio of the normal over in-plane magnitude 
of the external fleld, makes this magnetic hole system a 
good candidate for various applications, as the manip- 
ulation of large molecules using the magnetic holes, for 



example proteins which would be flxed to one or sev- 
eral holes coated with antigens, the determination of the 
transport properties of a ferrofluid, or as analog model 
of systems implying discrete particles of tunable inter- 
actions and hydrodynamic coupling to the carrier fluid. 
The effective pair interactions derived here should be a 
fundamental brick of any such applications. 
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